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We review the physics of charged impurities in the vicinity of graphene. The long-range nature of 
Coulomb impurities affects both the nature of the ground state density profile as well as graphene's 
transport properties. We discuss the screening of a single Coulomb impurity and the ensemble aver- 
aged density profile of graphene in the presence of many randomly distributed impurities. Finally, 
we discuss graphene's transport properties due to scattering off charged impurities both at low and 
high carrier density. 
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I. INTRODUCTION 

Graphene is a two dimensional sheet of carbon whose 
atoms arrange in a honeycomb lattice with nearest neigh- 
bor atoms forming strong sp2 bonds. The electronic 
properties of this material are mostly determined by the 
Pz orbitals with each carbon atom contributing one elec- 
tron to a Bloch band whose low energy properties are ade- 
quately described by a Dirac-Weyl effective Hamiltonian. 
While the study of Dirac Fermions has emerged in several 
contexts in theoretical condensed matter physics, its ex- 
perimental realization about three years ago, in the form 
of gated graphene devices where the carrier density 
can be tuned continuously from electron-like carriers for 
positive bias to hole-like carriers for negative gate volt- 
age, has prompted a prolific theoretical and experimental 
effort to understand the properties of this novel material. 

Most of excitement surrounding graphene stems from 
one of the following peculiar properties: (i) Electrons and 
holes in graphene have a gapless linear dispersion relation 
in contrast to the parabolic dispersion of other more con- 
ventional electron gases; (ii) The carriers in graphene are 
chiral - a property that has striking consequences such 
as the "half-integer" quantum Hall Effect^'^; and (iii) 
Carriers in graphene live at an exposed almost perfect 
2D surface that is amenable to surface probes^i^i^i^ and 
surface manipulation,^'^*' In addition, we note that there 
is the potential of mass producing graphene through epi- 
taxial growth methods^ and that graphene has remark- 
able mechanical propertiesi^iii^ which only further en- 
hance the interest. 

In this Perspective, we look at one important aspect 
of graphene which is the influence of disorder on its 
ground state and transport properties. We demonstrate 
that for graphene, charged (i.e. Coulomb) impurities be- 
have qualitatively different from neutral impuritiesi^ii^ii^ 
and dominate graphene's transport properties at low car- 
rier density. The importance of the Coulomb nature 
of graphene impurities was first highlighted by Ando^ii 
where by calculating the intraband contribution to the 
polarizability and absorbing the interband (i.e. electron- 
hole) contribution into a redefinition of the dielectric 
constant^ he showed that charged impurities could ex- 



plain the conductivity being linear-in-density as was seen 
in experiments^^ Similar conclusions were obtained by 
Nomura and MacDonald^° using a "complete screening" 
model (i.e. — > cxi), Cheianov and Falko using a nu- 
merical Thomas-Fermi approximation^^ and in Ref. [l^ 
using the full Random-Phase- Approximation (RPA) . An- 
alytic expressions for the RPA polarizability function cal- 
culated first in Ref. [H and then in Refs. [HlHIii re- 
vealed that for momentum transferred on the Fermi cir- 
cle (i.e. q = |k — k'l ~ 2kF sm9/2 < 2/cf) the graphene 
dielectric function calculated using the RPA was identi- 
cal to the much simpler Thomas-Fermi approximation 
at T = (see Fig. [1]). This then made it possible 
to calculate the RPA-Boltzmann conductivity analyti- 
cally,^^ and the dependence of graphene's conductivity on 
the fine-structure constant Ts = /{hvpi^) was recently 
verified experimentally.^" The importance of Coulomb 
scattering in explaining the observed graphene trans- 
port properties soon prompted an interest in investigat- 
ing the properties of a single charged impurity embedded 
in graphene. Katsnelson^ studied this problem using 
a Fermi- Thomas appro ximation, followed by studies in 

Refs. imOiiliSIlallllli who were mostly interested in 

effects beyond the RPA such as determining the critical 
impurity charge for which the Coulomb impurity forms 
bound states and the screening properties of graphene in 
the supercritical regime. 

It was understood by Refs. [2^13^ that as one ap- 
proached the Dirac point, one would soon encounter a 
situation where the gate voltage induced carrier density 
would be smaller than the fluctuation of carrier den- 
sity induced by the charged impurities thereby break- 
ing the graphene landscape into puddles of electrons and 
holes. Solving numerically for the conductivity using a 
finite-sized Kubo formalism for a limited range of impu- 
rity concentrations, Ref. [H concluded that the Coulomb 
disorder model gave a universal minimum conductivity 
whose value did not depend on the charged impurity con- 
centration, but that was larger than that expected for 
clean Dirac Fermions rMLi^^i^iiS while Ref. 12^ argued 
that this would give rise to a non-universal minimum 
conductivity whose value depended on the concentration 
of charged impurities. Ref. [13 developed a mean field 



2 



10 

8 
6 
4 
2 



2.5 

OJ 

1.5 


PDA 


^^^"^^ TF 

cs 


I 


1 1.5 

X 







1 2 3 4 5 6 7 8 
x= q/2kp 

FIG. 1: The main panel shows different dielectric functions 
used in the literature, including the "Complete Screening" 
(CS), "Thomas- Fermi" (TF) and "Random Phase Approxima- 
tion" (RPA). The inset shows a blow-up at q = 2fcF to show 
where "Step Approximation" (SA) used in Ref.[23 differs from 
the exact result. 



approach to understand the properties of graphene at 
the Dirac point by calculating an effective carrier den- 
sity self-consistently. This theory made quantitative pre- 
dictions about the dependence of the minimum conduc- 
tivity and rms carrier density on the charged impurity 
concentration and substrate dielectric constant, and in 
particular argued that cleaner graphene samples would 
have larger minimum conductivity. Ref. ^47 then stud- 
ied the ground state properties of graphene by mini- 
mizing an energy functional comprising kinetic energy, 
Hartree, exchang o^^'^^ i^^ and corrclation^^i^i^ contribu- 
tions in the presence of Coulomb disorder. This work 
made quantitative predictions about properties of the 
carrier density distribution, both at and away from the 
Dirac point, and enabled Ref. [4^ to develop an effec- 
tive medium theory to calculate the graphene's conduc- 
tivity through these inhomogeneous puddles, capturing 
quantitatively the minimum conductivity plateau that is 
seen in experiments^^ '^^i^^ We mention that underlying 
the existence of this minimum conductivity plateau is 
the high transmission of graphene p-n junctions, which 
has been the subject of theoretical^i^ and experimental 
studyi^i^i^ For the purposes of this paper we do not 
discuss quantum interference effects (see Ref. [s^ and ref- 
erences therein) or the strongly interacting regime (see 
Ref. m and references therein) . 

The remainder of this paper is structured as follows. 
In Section [IT] we discuss the problem of the screening 
of a single Coulomb impurity in the sub-critical regime 
as a useful toy model to understand the many impurity 
problem that we address in Section IIIII where we study 
the case of many Coulomb impurities that are uncor- 
related and distributed uniformly in order to study the 
ground state properties of graphene. In Section IIV Al 



we review the high-density Boltzmann transport theory, 
and discuss the Effective Medium Theory (EMT) in Sec- 
tion IIVBI In Sections |Vl IVIi and IVIII we briefly review 
the experimental situation, discuss graphene minimum 
conductivity, and recent theoretical work not covered in 
this review. We then conclude in Section IVlIII 



II. SCREENING OF A SINGLE COULOMB 
IMPURITY 

Following Ref. [H one can construct the Thomas-Fermi 
screening of a single charged impurity. The goal is to 
calculate the screened Coulomb potential Vs{r) = Vo{r) + 
Vlnd(''), where the bare potential Vq = hv-prs[r'^ + 
and the induced potential is given by 
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Vind{r) = (hvprs) / dr' 



(1) 



where we can imagine tuning the back gate to ensure 
charge neutrality n ~ Q. If one further assumes that 
the local carrier density is given by the Fermi-Thomas 
condition, n[T^(r)] = —V{rY/['K{hvF)'^], one can write 
down a (one dimensional) self-consistency equation for 
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where K[x\ is the complete elliptic integral of the first 
kind. The screened potential induced for this single im- 
purity using this method was discussed in Ref. |2a. This 
formalism can be generalized using the method devel- 
oped in Ref. |4ll to include the effects of exchange. The 
ground state carrier density can be obtained from the 
Thomas-Fermi-Dirac (TFD) energy functional 
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where the first term in Eq. [3] is the kinetic energy, the 
second term is the Hartree part of the Coulomb interac- 
tion, the third is the exchange-correlation energy and the 
fourth term is the energy due to disorder, where Vd is 
the disorder potential and the last term is added to set 
the average carrier density, (n) , through the chemical po- 
tential A. The correlation term is much smaller than ex- 
change and, to very good approximation^^*^!^ is propor- 
tional to exchange. Therefore, hereafter, we neglect the 
correlation contribution by assuming SExc/Sn = I](n), 
where S(rt) is the Hartree-Fock self-energy.— i^i^ The 
energy functional Eq. [3] is quite general and can be tai- 
lored by properly choosing Vd and its coupling to n{r), 
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to consider different sources of disorder. For the single 
impurity problem, the solution can also be cast as a one 
dimensional integral equation 
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where A = 1/(0.25 nm) is the band energy cutoff and 
C « 0.916. 



III. GROUND STATE PROPERTIES AT THE 
DIRAC POINT 

The single impurity problem discussed in the previ- 
ous section is a much simpler problem because rotational 
symmetry makes the problem one-dimensional. Adding 
many impurities also brings further complications: while 
the carrier density induced by a single impurity is neg- 
ligible, this is not the case for many impurities where 
although the average density can be tuned to zero via an 
external gate potential, the scale of the density fluctua- 
tions is set by the impurity concentrationj^ As shown in 
Fig. [H the RPA screening properties of graphene are very 
different at the Dirac point (i.e. fcp — ^ 0) and at finite 
density, therefore theoretical frameworks constructed to 
work at the Dirac point are bound to fail when there 
are such large density fluctuations. In this section, we 
present two different approaches to describe the Dirac 
point. The first is a mean- field theory where an effective 
density n* is obtained by solving self-consistently for the 
density induced by the fluctuations of the screened impu- 
rity potential (that itself depends on density). The sec- 
ond is a generalization of the energy functional method 
discussed above for a single impurity to the much more 
complicated case of many Coulomb impurities. 



A. Self-consistent Approximation (SCA) 




FIG. 2: Spatial correlation function for the screened poten- 
tial at the Dirac point for Vs = 0.8 and d = 1 nm. Red 
diamonds are the results obtained by minimizing graphene's 
energy functional (Sec. IIIIB|I . The lines are the SCA re- 
sults using Eq. I5al blue lines, and its Gaussian approxima- 
tion, Eq. I5bl All the results are normalized via the value of 
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Fermi energy {V'^)c = {Epln*])'^ = TT{hvpYn* . This self- 
consistent approximation then allows us to compute any 
correlation function at the Dirac point, although closed 
form analytic results are often elusive. To make analyti- 
cal progress, one can map 
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where analytic expressions for Kq and ^, were reported in 
Ref. '56^. A numerical evaluation of Eq. [5a| and the Gaus- 
sian approximation Eq. I5bl is shown in Fig. [51 Within 
the Gaussian approximation one finds that nj-ms = 
^ly^/[TTihvF)^] w n* ^3 + [mmpne]-' « V3n*, where 
in the last equation we further assume that ni,„p7r^^ ~ 
^ 1. This result for rij-^s is particularly useful when 
comparing the self-consistent approximation with other 
methods. 



For any microscopic single impurity potential (f){r,n), 
the probability distribution, P{V), of the total poten- 
tial, V , is P{V) = {5{V ~ J2flT Hn,ri)))r^ where 
{■ ■ ■)ri is the average over all possible disorder configura- 
tions. Assuming that the impurities positions are uncor- 
related one can compute expressions for all moments of 
the induced disorder potentiali^i^ For example, the con- 
nected moment {V'"}c = riimp / d^r[(/)[(r, n)]*-'. The self- 
consistent approximation involves obtaining the effective 
carrier density n* by equating the second moment of the 
disorder potential with the square of the corresponding 



B. Energy Functional Minimization (EFM) 

To study graphene transport properties for a distribu- 
tion of charged impurities, we use Eq. [3] taking Vd to 
be the potential generated by a random 2D distribution 
C(r) of impurity charges placed at a distance d from 
the graphene layer. We assume C(r) to be on average 
zero and uncorrelated and perform our calculations on 
a 200 nmx200 nm square sample with a 1 nm spatial 
discretization. Close to the Dirac point, for a single dis- 
order realization, we find that the carrier density breaks 



4 



up into electron-hole puddles. Since we are interested in 
disorder averaged quantities, we examine several disorder 
realizations (500-1000) and denote disordered averaged 
quantities by angled brackets. To characterize the den- 
sity profile, we calculate the disorder averaged density- 
density correlation fimction {dn{r)6n{0)) , from which we 
can extract the root mean square rirms = \/ {Sn{0)Sn{0)) , 
and the typical correlation length, ^, defined in this sec- 
tion as the FWHM of {Sn{r)5n{0)) . We find^^ for typ- 
ical graphene samples, that rirms ~ (n) for dopings as 
high as 10^^ cm"^ and that close to the Dirac point. 



for 



''imp 



< 10^" cm- 



including exchange is three 



times smaller than without. We find the correlation 
length ^ to be of the order of 10 nm, see Fig. [31 This 
value suggest that the electron-hole puddles are quite 
small. However a closer inspection reveals that, close 
to the Dirac point, the density profile is characterized by 
two distinct types of inhomogeneities^'^: wide regions (i.e. 
big puddles spanning the system size) of low density con- 
taining a number of electrons (holes) of order 10; and few 
narrow regions, whose size is correctly estimated by ^, of 
high density containing a number of carriers of order 2. 
This picture is confirmed by the results shown in Fig. [4] in 
which the disorder averaged area fraction, Aq, over which 
|n(r) — < rirms/ W is plotted as a function of riimp- 
We see that is of order 1 /3 and we also find that the 
area fraction over which |n(r) — (n)| is less than 1/5 of 



is close to 50% for n. 



imp 



The 



bination of the relatively high density in the peaks/dips 
and the fact that in the low density regions n(r) varies 
over scales much bigger than 10 nm guarantees that the 
inequality ■y/7rn[| Vn|/n]~"^ ^ 1 is satisfied over the ma- 
jority of the graphene sample and therefore justifies the 
use of the EFM theory. The EFM should be a reasonable 
quantitative theory for existing graphene samples at all 
values of the carrier density. 
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FIG. 3; Density correlation length ^ at the Dirac point as 
function of nimp for Vs = 0.8 and two diflerent values of d. 
Results with (without) exchange are shown by solid (dashed) 
lines. 
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C. Comparison of SCA and EFM 

Here we compare the results from the Self- Consistent 
Approximation of Sec. IIII Al and the Energy Functional 
Minimization (EFM) formalism discussed in Sec IIIIBI 
Figure [2] shows the disordered averaged spatial correla- 
tion function at the Dirac point for the screened disorder 
potential V = Vd + (1/2) f (fr'n{r')/\r - r'\ where the 
(red) diamonds are obtained minimizing Eq. [3l The solid 
(blue) line shows the same quantity calculated using the 
self-consistent approximation (SCA). The EFM approach 
and the SCA give a similar behavior for {V(r)V{0)), char- 
acterized by an algebraic, cx decay at large dis- 
tances. The green solid line shows a Gaussian approxi- 
mation which captures much of the quantitative details 
of the screened disorder potential correlation function, 
but not the power law 1/r^ decay. 

In Figs. 15171 rirms at the Dirac point is shown as func- 
tion of riimp , Ts and d respectively. The red dashed lines 
show the results obtained using the EFM theory includ- 



FIG. 4: Area fraction, Aq, over which is |n(r) — {n)\ < 
rirms /lO as a function of at the Dirac point for Vs = 0.8. 

Results with (without) exchange are shown by solid (dashed) 
lines. 

ing exchange, the solid blue lines are the results obtained 
using the SCA theory. In general the SCA gives values of 
rirms smaller than the EFM theory but in general there 
is good semi-quantitative agreement especially at low rg 
and ri-irrip • 

IV. GRAPHENE CONDUCTIVITY 

A. High-density: Boltzmann Transport Theory 

In this section we investigate the graphene transport 
for large carrier densities (n ^ n^), where the system is 
homogeneous. We show in detail the microscopic trans- 
port properties at high carrier density using the Boltz- 
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FIG. 5: Results for rirms as a function of Wimp at the Dirac 
point for rs — 0.8. The red dashed hues are the results ob- 
tained minimizing graphene's energy functional and the solid 
lines are the SCA results. 
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FIG. 7: Resuhs for function of d at the Dirac point 

for riimp = 10^^ cm~^ and rs — 0.8. The red dashed lines are 
the results obtained minimizing graphene's energy functional 
and the solid lines are the SCA results. 
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FIG. 6: Resuhs for function of rs at the Dirac point 

for two values of nimp and d = 1 nm. The red dashed lines are 
the results obtained minimizing graphene's energy functional 
and the solid lines are the SCA results. 



mann transport theoryi^ We calculate the conductivity 
a (or mobility fi = a/ne) in the presence of randomly 
distributed Coulomb impurity charges near the surface 
with the electron-impurity interaction being screened by 
the 2D electron gas in the random phase approximation 
(RPA). Even though the screened Coulomb scattering is 
the most important scattering mechanism in our calcu- 
lation, there are additional scattering mechanisms (i.e. 
neutral point defects) unrelated to the charged impu- 
rity scattering for very high mobility samples. Point de- 
fects gives rise to a constant conductivity in contrast to 
charged impurity scattering which produces a conductiv- 
ity linear in n/rii. Our formalism can include both effects, 
where zero range scatterers are treated with an effective 
point defect density of Up. For the purpose of this cal- 
culation, we neglect all phonon scattering effects, which 



were considered recently with the finding that acoustic 
phonon scattering gives rise to a resistivity that is linear 
in temperatureJ^S 

We start by assuming graphene to be a homogeneous 
2D carrier system of electrons (or holes) with a carrier 
density n induced by the external gate voltage. The low- 
energy band Hamiltonian for homogeneous graphene is 
well-approximated by a 2D Dirac equation for massless 
particles, 

H = hvF{(Txkx + <7yky), (6) 

where vp the 2D Fermi velocity, ax and ay are Pauli 
spinors and k is the momentum relative to the Dirac 
points. The corresponding eigenstates are given by the 
plane wave Vsk(r) = ^exp(ik • r)Fsu, where A is 
the area of the system, s = ±1 indicate the conduc- 
tion (-fl) and valence (—1) bands, respectively, and 

^Jk ~ "Tf (^'^''' '^^ with 6'k — ta.n{ky/kx) being the po- 
lar angle of the momentum ftk. The corresponding 
energy of graphene for 2D wave vector k is given by 
Esk = sfi.WF|k|, and the density of states (DOS) is given 
by D{e) = g\e\/{2TTh'^Vp), where g = QsQv is the total 
degeneracy {gs = 2,(7„ = 2 being the spin and valley 
degeneracies, respectively). 

When the external force is weak and the displacement 
of the distribution function from the thermal equilibrium 
value is small, we can use linearized Boltzmann equation 
within relaxation time approximation. In this case the 
conductivity for graphene can be written by 

o- = j dekD{ek)T{ek) ( I . (7) 

Note that /(cfc) is the Fermi distribution function, 
f{ek) = {1 + exp[(efc - \)]/kBT}-^ where the finite 
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temperature chemical potential A(T) is determined self- 
consistently to conserve the total number of electrons. 
At T = 0, /(e) is a step function at the Fermi energy 
Ep = X{T ~ 0), and we then recover the Einstein re- 
lation a = ^D{Ef)t{Ef). In Eq. © T(esk) is the 
relaxation time or the transport scattering time of the 
collision and is given by 
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ft ^'"^ J (27r)2l^^k,.k' 
[1 - cos^kk-lt^ (csk - e^k') , 



(8) 



where 0kk' is the scattering angle between the scattering 
in- and out- wave vectors k and k', and {V^■^\,-^^,) is the 
matrix element of the scattering potential associated with 
impurity disorder in the graphene environment, and n^"'' 
is the number of impurities per unit area of the a-th 
kind of impurity. Note that since we consider the elastic 
impurity scattering the interband processes (s ^ s') are 
not permitted. 

The matrix element of the scattering potential of ran- 
domly distributed screened impurity charge centers in 
graphene is given by 



eiq) 



l-f cos 6* 
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where q — |k — k'|, 6 = 6'kk', and Vi{q,d) = 
27Te^ exp{—qd) / (Kq) is the Fourier transform of the 2D 
Coulomb potential in an effective background lattice di- 
electric constant k, where d is the location of the charged 
impurity measured from the graphene sheet. The factor 
(l-fcos0)/2 arises from the sublattice symmetry (overlap 
of wave function) fii, 

In Eq. (O, e{q) = e{q,T) is the 2D finite temper- 
ature static RPA dielectric (screening) function appro- 
priate for graphene, given by e{q,T) = 1-1- Vc{q)[^ — 
G{q)]II{q,T), where U{q,T) is the graphene irreducible 
finite-temperature polarizability function;^ Vc{q) is the 
Coulomb interaction, and G{q) is the local field correc- 
tion. In RPA, G{q) = and in Hubbard approximation 
(HA), G{q) 



In Fig. [8] we show the calculated graphene conductiv- 
ities limited by screened charged impurities. The RPA 
screening used in our calculation is the main approxi- 
mation. We also show results for HA screening which 
includes local field corrections approximately. We note 
that since the graphene is a weakly interacting system 
(rs < 1) the correlation effects are not strong. We em- 
phasize that in order to get quantitative agreement with 
experiment, the screening effects must be included. Using 
the unscreened dielectric function would have conductiv- 
ity less than Ae^/h for the entire range of gate voltages 
used in the experiment. Our main result with screened 
Coulomb impurities is the quantitative agreement with 
experiments in the regime where the conductivity is lin- 
ear in density. In inset we show the effect of remote 
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FIG. 8: Calculated graphene conductivity as a function of 
carrier density {n/rn, where rn is the impurity density) lim- 
ited by Coulomb scattering. RPA (HA) indicates the result 
with the RPA (HA) screened Coulomb scattering, and "no 
screening" indicates the results with bare Coulomb scatter- 
ing. Note that the calculated conductivity with unscreened 
Coulomb potential is less than 4e^//i for given density range. 
In inset the effect of remote scatters is shown. Here d is the 
distance between the 2D graphene layer and the 2D impurity 
layer. 



scatterers which are located at a distance d from the in- 
terface. The main effect of remote impurity scatterings 
is that the conductivity deviates from the linear behavior 
with density and increases with both the distance d and 
n/n.i. 

For very high mobility samples, one finds a sub-linear 
conductivity instead of the linear behavior with den- 
sity. Such high quality samples presumably have a small 
charge impurity concentration Ui and it is therefore likely 
that point defects here play a more dominant role. Point 
defects gives rise to a constant conductivity in contrast to 
charged impurity scattering which produces a conductiv- 
ity linear in n/n^ . In the presence of both the long-ranged 
charged impurity and the short-ranged neutral impurity, 
the total scattering time becomes l/n ~ 1/rj + I/tq 
where (tq) is the scattering time due to charged 
Coulomb (short ranged) impurities. Shown in Fig. [9] 
is the graphene conductivity calculated including both 
charge impurity and zero range point defect scattering 
for different ratios of the point scatterer impurity density 
Up and the charge impurity density n;. For small np/ui 
we find the linear conductivity that is seen in most ex- 
periments and for large Up/rii we see the flattening out of 
the conductivity curve (which in the literature^ has been 
referred to as the sub- linear conductivity). We believe 
this high-density flattening of the graphene conductiv- 
ity is a non-universal crossover behavior arising from the 
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n/ri: 



FIG. 9: Graphene conductivity calculated using a combina- 
tion of short and long range scatterers. In this calculation, 
we use np/rii = 0, 0.01, 0.02 (top to bottom). In inset we 
show the graphene mobility as a function of dielectric con- 
stant (k) of substrate for different carrier densities n = 0.1, 
1, 5 X lO^^cm"^ (from top to bottom) in the presence of both 
long ranged charged impurity {rii = 2 x lO^'^ cm~^) and short- 
ranged neutral impurity {rip — 0.4 x 10^" cm~^). Vb = 1 
KeVA^ is used in this calculation, which corresponds to the 
Coulomb potential of electron density n — 10^^ cm~^. 
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FIG. 10: Calculated conductivity for different temperatures 
T = 0, 200, 300, 500K (top to bottom) as a function of density 
with an impurity density = 5 x lO^^cm"^. 



competition between two kinds of scatterers. In general 
this crossover occurs when two scattering potentials are 
equivalent, that is, riiV^ ~ npV^ . 

In the inset of Fig. [H] we show our calculated mobility 



in the presence of both charged impurities and short- 
ranged impurities. As the scattering limited by the short- 
ranged impurity dominates over that by the long-ranged 
impurity (e.g. npVi^ 3> n-iV^) the mobility is no longer 
dependent on the charged impurity and approaches its 
limiting value 
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The limiting mobility depends only on neutral impurity 
concentration Up and carrier density, which indicates that 
to get high graphene mobility it is necessary to have de- 
fect free graphene. 

Finally in Fig. [10] we show the calculated temperature 
dependent conductivity for different temperatures as a 
function of density. We note that there are two inde- 
pendent sources of temperature dependent resistivity in 
our calculation. One comes from the energy averaging 
defined in Eq. ([7]), and the other is the expHcit temper- 
ature dependence of the dielectric function e{q, T) which 
produces a temperature dependent t{s,T). Figure [TO] 
shows that in the high density limit the conductivity 
decreases as the temperature increases, but in the low 
density limit the conductivity shows non-monotonic be- 
havior, i.e. cr(T) has a local minimum at a finite temper- 
ature and increases as the temperature increases. Thus, 
we find that the calculated conductivity shows a non- 
monotonicity in the low density limit, i.e., at low temper- 
atures the conductivity shows metallic behavior and at 
high temperatures it shows insulating behavior. The non- 
monotonicity of temperature dependent cr(T) is under- 
stood to arise from temperature dependent screeningi^ 
We mention that for T > 100 K, phonons contribute 
to the temperature dependence of graphene conductiv- 



B. Low Density: Effective Medium Theory 

At low density the fluctuations in carrier density be- 
come larger than the average density. To understand 
the transport properties of this inhomogeneous system, 
Ref. developed an effective medium theory where 
graphene's conductivity is found by solving an integral 
equation 



a[n\ + CTEMT 



(11) 



where P[n\ is the density distribution function and 
a[n\ is the (local) Boltzmann conductivity discussed in 
Sec. IIV Al For the purpose of this section, we take 
o'bM = (2e^/G[r's]/i)n/nii„p where G[rs] was derived in 
Ref.[27land shown in Fig. [11] The quantitatively accurate 
theory using P[n] derived from the EFM of Sec lIII B] was 
developed in Ref. 45. Here we derive analytical results 
obtained by using model distribution functions for P[n], 
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FIG. 11: The main panel shows G[rs] that parameterizes the 
inverse scattering time in the Boltzmann theory. The analytic 
form of G[rs] can be found in Ref. [23- The dashed line shows 
the experimentally relevant regime for graphene on Si02 sub- 
strates. Inset shows sample mobility as a function of substrate 
dielectric constant Ks for nimp — 2x 10^^ cm~^. Changing Ks 
by a factor of 2, increases mobility by 50 percent. 



which as discussed in Ref. 45 show quantitative agree- 
ment with the numerical theory only for smaU and low 
nimp- To illustrate this method, we first consider P[n\ to 
be a Gaussian distribution. Requiring that J n'^P[n] = 
n^uis fixes all the free parameters. Solving Eq. [11] then 
gives zexp~^'(7rErfi[z] — Ei[z^]) — y/TT/2, where Erfi is 
the imaginary error function, Ei is the exponential inte- 
gral function and z = crEMT/(\/2a-_B [nrms]) ~ 0.405 giv- 
ing (Temt ~ 0.9925 ctscAj where we use the results of 
Sec. lIII Al that (Tsca = o'B[rt*] and rtj-ms ~ V^n* . The de- 
velopment of an effective medium theory for graphene^ 
now allows us to reinterpret the results of Ref. [l^ as 
equivalent to the assumption that P[n] is Gaussian with 
density fluctuations determined by the self-consistency 
condition = (Vd)^. 

We can explore other functional forms for P[n]. For 
a Lorentzian P[n] = (nz,/7r)/(n| -I- n^), one can solve 
Eq. [11] analytically giving (Temt — crsinL]- If one 
identifies the width of the Lorentzian with the self- 
consistent carrier density til = n* , then this provides 
another way to understand the self-consistent transport 
result. Subsequent to the results of Ref. ^451, Fogler 
developed an effective medium theory using P[n] — 

(1/V2 ) exp[-\/2|n|/nrms], where similar to the 
Gaussian distribution case, requiring normalization and 
setting {n?) — ri^^^^ fixes all the free parameters. Solv- 
ing Eg. [TT] gives ze^r[z] = 1/2 where F is the Gamma 
function and z = \/2crEMT/o'_B[nrms] ~ 0.610. Again, ap- 
proximating firms ~ -s/Sn*, we find ctemt ~ 0.75 ctsca, 
which is different from the numerical results obtained in 
Ref. ig^. While these analytical approximations are use- 
ful in providing a qualitative understanding of graphene 
transport, quantitative differences remain between these 
and the full numerical solution'*^ especially at large im- 



purity concentrations and large (See e.g. Fig. [6]). 



C. Suspended Graphene 

One of the direct consequences of charged impurity 
scattering in graphene is the prediction32i that the elimi- 
nation of charged impurities from the graphene environ- 
ment, for example, by suspending graphene, without any 
substrate would lead to a much enhanced carrier mobility. 
Recently, Bolotin et a L^^'^'^ managed to remove charged 
impurities from the graphene environment by suspend- 
ing graphene without any substrate (and simply current- 
annealing away any remnant impurities on the graphene 
surface). This immediately led to an oder of magnitude 
increase (to /i ^ 10^ cm^/Vs) in the graphene mobil- 
ity as predicted theoretically. Recent theoretical work^ 
shows excellent agreement with the transport measure- 
ments on suspended graphene^^i^i^ with both the re- 
duced impurity density and the modified screening (due 
to the elimination of the substrate) contributing to the 
graphene conductivity. Since phonon scattering effects 
in graphene are weak upto room temperature,^- the en- 
hanced graphene mobility arising from the elimination of 
charged impurities may lead to very high (> 10^ cm^/Vs) 
graphene mobilities even at room temperature. 



V. DISCUSSION OF EXPERIMENTS 

One of the very first puzzles in graphene transport 
experiments^ was that the conductivity was linear in 
carrier density, whereas existing theory'^^ predicted con- 
stant conductivity at high density. As discussed in 
the introduction, the linear in density emerges naturally 
from the Boltzmann transport theory of charged impu- 
rities j^i^i^'^^i^l and to our knowledge, no other the- 
ory produces this linear behavior without a fine-tuning 
of parameters to make the scattering mimic Coulomb 
impurities. Moreover, three recent experiments have 
rigorously verified the high density predictions of the 
Boltzmann transport theory. In Ref. [l9| sample mobility 
was correlated with shift of the Dirac point and plateau 
width showing qualitative and semi-quantitative agree- 
ment with the theory presented here. Reference l46l were 
able to directly measure the effect of Coulomb scatterers 
by intentionally adding potassium ions to graphene in 
ultra-high vacuum observing qualitatively all the predic- 
tions of the (self-consistent) Boltzmann theory. Finally, 
Ref. [13 was able to tune graphene's fine structure con- 
stant by depositing ice on top of graphene, and thereby 
precisely testing predictions of the Boltzmann theory. 
Results from this experiment are shown in Fig 1121 Since 
the dielectric constant of the Si02 substrate and that of 
ice are known, there are no adjustable parameters in the 
theoretical curve. 
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FIG. 12: Data taken from Ref. |lO| shows the effect of di- 
electric screening on graphene. Data points show the dif- 
ference in graphene conductivity before and after depositing 
ice, while solid line shows the theory with no adjustable pa- 
rameters. The non-monotonic behavior is a consequence of 
the competing effects of dielectric screening on Coulomb and 
short-range scatterers, whereas the quantitative agreement is 
a stringent test of the the Boltzmann theory. Inset, also taken 
from Ref. [T^ shows the raw experimental data. 



the Boltzmann transport theory provides excellent de- 
scription of the non-universal behavior of the minimum 
conductivity around the Dirac point. In particular, the 
sample dependence of the minimum conductivity arises 
from the different impurity disorder in different samples. 

Quantum effective field theories of graphene minimum 
conductivity, which predict a universal minimum conduc- 
tivity, are inapplicable to real graphene samples because 
real graphene is dominated by disorder induced inhomo- 
geneity near the Dirac point, which is outside the scope of 
the quantum field theories. Although the self-consistent 
effective medium theory developed by us^lii^ gives rea- 
sonable agreement with the experimental observations, 
the key conceptual question of what happens at T = as 
disorder also goes to zero still remains open. Such a sce- 
nario is of course experimentally irrelevant (since experi- 
ments are performed at finite temperatures in disordered 
systems), but the theoretical question of the conductiv- 
ity crossover from the inhomogeneous Boltzmann regime 
of Refs. [27ll45l to the homogeneous quantum transport 
regime is an interesting open question. Recent attempts 
to understand such crossover phenomena include several 
complementary theoretical avenues (See Refs. [Ulle^ 
and references therein). 



VII. RECENT WORK 



VI. THE MINIMUM CONDUCTIVITY PUZZLE 

As described above in Sec. IIVBI our recent theoreti- 
cal work«^2iM provides a satisfactory explanation for the 
minimum conductivity phenomenon in graphene near the 
charge neutrality (i.e. Dirac) point. In particular, early 
theoretical work^Sii^ii^i^SdS predicted a universal T = 
minimum conductivity CTmin — '^e^ /irh at the graphene 
Dirac point in clean disorder-free systems. The inclusion 
of disorder- induced quantum anti- localization effect, as- 
suming no intervalley scattering, leads to a theoretical in- 
finite minimum conductivity at the Dirac point, whereas 
the presence of inter-valley scattering localizes the sys- 
tem leading to zero conductivity at the Dirac point. This 
confusing theoretical picture stands in stark contrast to 
the experimental reality, where the graphene conductiv- 
ity is approximately a constant (as a function of gate volt- 
age or carrier density) around the Dirac point, with this 
constant minimum conductivity plateau having a non- 
universal sample dependent value (~ ie^/h — 20e^//i). 

It was first suggested in Ref. [l^ that the minimum con- 
ductivity phenomenon is closely related to the break-up 
of the graphene landscape into inhomogeneous puddles 
of electrons and holes around the Dirac point due to the 
effect of the charged impurities in the environment. This 
physical idea was further developed into a quantitatively 
successful theory (see Sec. IIVBI above) in Refs. [27ll45l . 
where it was shown that a self-consistent treatment of 
the impurity induced electron-hole puddles coupled with 



Among recent relevant work not discussed in this re- 
view we mention a detailed calculation of the temper- 
ature dependent graphene conductivity due to electron- 
phonon scattering;^ a detailed calculation of the temper- 
ature dependent graphene conductivity due to the tem- 
perature dependence of the screening of charged impurity 
scattering^ a calculation of graphene density of states 
as modified by impurity scattering,®^ a consideration of 
percolation induced localization transition in graphene 
nanoribbons,^® a theory of charged impurity screening in 
graphene bilayers,®- and a prediction of graphene magne- 
toresistance induced by a parallel magnetic field through 
the spin-polarization dependence of screeningi^ 



VIII. CONCLUSION 

We have developed here the theory for Coulomb im- 
purities on graphene. As we have shown. Coulomb im- 
purities behave qualitatively different from short-range 
scatterers such as point defects or missing atoms. Away 
from the Dirac point, the physics is well described by a 
semi-classical Boltzmann transport theory, while at the 
Dirac point density fluctuations dominate breaking the 
system into puddles of electrons and holes. We have 
shown how these inhomogeneities can be characterized 
by a mean-field self-consistent theory and by numerically 
minimizing graphene's energy functional, and that an ef- 
fective medium theory can be employed to describe the 
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low density transport properties giving semi-quantitative 
agreement with experimental results. 
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